home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / claqhe.z / claqhe
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. CCCCLLLLAAAAQQQQHHHHEEEE((((3333FFFF))))                                                          CCCCLLLLAAAAQQQQHHHHEEEE((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CLAQHE - equilibrate a Hermitian matrix A using the scaling factors in
  10.      the vector S
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE CLAQHE( UPLO, N, A, LDA, S, SCOND, AMAX, EQUED )
  14.  
  15.          CHARACTER      EQUED, UPLO
  16.  
  17.          INTEGER        LDA, N
  18.  
  19.          REAL           AMAX, SCOND
  20.  
  21.          REAL           S( * )
  22.  
  23.          COMPLEX        A( LDA, * )
  24.  
  25. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  26.      CLAQHE equilibrates a Hermitian matrix A using the scaling factors in the
  27.      vector S.
  28.  
  29.  
  30. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  31.      UPLO    (input) CHARACTER*1
  32.              Specifies whether the upper or lower triangular part of the
  33.              Hermitian matrix A is stored.  = 'U':  Upper triangular
  34.              = 'L':  Lower triangular
  35.  
  36.      N       (input) INTEGER
  37.              The order of the matrix A.  N >= 0.
  38.  
  39.      A       (input/output) COMPLEX array, dimension (LDA,N)
  40.              On entry, the Hermitian matrix A.  If UPLO = 'U', the leading n
  41.              by n upper triangular part of A contains the upper triangular
  42.              part of the matrix A, and the strictly lower triangular part of A
  43.              is not referenced.  If UPLO = 'L', the leading n by n lower
  44.              triangular part of A contains the lower triangular part of the
  45.              matrix A, and the strictly upper triangular part of A is not
  46.              referenced.
  47.  
  48.              On exit, if EQUED = 'Y', the equilibrated matrix:  diag(S) * A *
  49.              diag(S).
  50.  
  51.      LDA     (input) INTEGER
  52.              The leading dimension of the array A.  LDA >= max(N,1).
  53.  
  54.      S       (input) REAL array, dimension (N)
  55.              The scale factors for A.
  56.  
  57.      SCOND   (input) REAL
  58.              Ratio of the smallest S(i) to the largest S(i).
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. CCCCLLLLAAAAQQQQHHHHEEEE((((3333FFFF))))                                                          CCCCLLLLAAAAQQQQHHHHEEEE((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      AMAX    (input) REAL
  75.              Absolute value of largest matrix entry.
  76.  
  77.      EQUED   (output) CHARACTER*1
  78.              Specifies whether or not equilibration was done.  = 'N':  No
  79.              equilibration.
  80.              = 'Y':  Equilibration was done, i.e., A has been replaced by
  81.              diag(S) * A * diag(S).
  82.  
  83. PPPPAAAARRRRAAAAMMMMEEEETTTTEEEERRRRSSSS
  84.      THRESH is a threshold value used to decide if scaling should be done
  85.      based on the ratio of the scaling factors.  If SCOND < THRESH, scaling is
  86.      done.
  87.  
  88.      LARGE and SMALL are threshold values used to decide if scaling should be
  89.      done based on the absolute size of the largest matrix element.  If AMAX >
  90.      LARGE or AMAX < SMALL, scaling is done.
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.